Improved production of the antidiabetic metabolite montbretin A in Nicotiana benthamiana: discovery, characterization, and use of Crocosmia shikimate shunt genes.

The plant-specialized metabolite montbretin A (MbA) is being developed as a new treatment option for type-2 diabetes, which is among the ten leading causes of premature death and disability worldwide. MbA is a complex acylated flavonoid glycoside produced in small amounts in below-ground organs of the perennial plant Montbretia (Crocosmia × crocosmiiflora). The lack of a scalable production system limits the development and potential application of MbA as a pharmaceutical or nutraceutical. Previous efforts to reconstruct montbretin biosynthesis in Nicotiana benthamiana (Nb) resulted in low yields of MbA and higher levels of montbretin B (MbB) and montbretin C (MbC). MbA, MbB, and MbC are nearly identical metabolites differing only in their acyl moieties, derived from caffeoyl-CoA, coumaroyl-CoA, and feruloyl-CoA, respectively. In contrast to MbA, MbB and MbC are not pharmaceutically active. To utilize the montbretia caffeoyl-CoA biosynthesis for improved MbA engineering in Nb, we cloned and characterized enzymes of the shikimate shunt of the general phenylpropanoid pathway, specifically hydroxycinnamoyl-CoA: shikimate hydroxycinnamoyl transferase (CcHCT), p-coumaroylshikimate 3'-hydroxylase (CcC3'H), and caffeoylshikimate esterase (CcCSE). Gene expression patterns suggest that CcCSE enables the predominant formation of MbA, relative to MbB and MbC, in montbretia. This observation is supported by results from in vitro characterization of CcCSE and reconstruction of the shikimate shunt in yeast. Using CcHCT together with montbretin biosynthetic genes in multigene constructs resulted in a 30-fold increase of MbA in Nb. This work advances our understanding of the phenylpropanoid pathway and features a critical step towards improved MbA production in bioengineered Nb.

Genomic virulence features of Beauveria bassiana as a biocontrol agent for the mountain pine beetle population.

BACKGROUND: The mountain pine beetle, Dendroctonus ponderosae, is an irruptive bark beetle that causes extensive mortality to many pine species within the forests of western North America. Driven by climate change and wildfire suppression, a recent mountain pine beetle (MPB) outbreak has spread across more than 18 million hectares, including areas to the east of the Rocky Mountains that comprise populations and species of pines not previously affected. Despite its impacts, there are few tactics available to control MPB populations. Beauveria bassiana is an entomopathogenic fungus used as a biological agent in agriculture and forestry and has potential as a management tactic for the mountain pine beetle population. This work investigates the phenotypic and genomic variation between B. bassiana strains to identify optimal strains against a specific insect. RESULTS: Using comparative genome and transcriptome analyses of eight B. bassiana isolates, we have identified the genetic basis of virulence, which includes oosporein production. Genes unique to the more virulent strains included functions in biosynthesis of mycotoxins, membrane transporters, and transcription factors. Significant differential expression of genes related to virulence, transmembrane transport, and stress response was identified between the different strains, as well as up to nine-fold upregulation of genes involved in the biosynthesis of oosporein. Differential correlation analysis revealed transcription factors that may be involved in regulating oosporein production. CONCLUSION: This study provides a foundation for the selection and/or engineering of the most effective strain of B. bassiana for the biological control of mountain pine beetle and other insect pests populations.

4-Coumaroyl-CoA ligases in the biosynthesis of the anti-diabetic metabolite montbretin A.

BACKGROUND: Montbretins are rare specialized metabolites found in montbretia (Crocosmia x crocosmiiflora) corms. Montbretin A (MbA) is of particular interest as a novel therapeutic for type-2 diabetes and obesity. There is no scalable production system for this complex acylated flavonol glycoside. MbA biosynthesis has been reconstructed in Nicotiana benthamiana using montbretia genes for the assembly of MbA from its various different building blocks. However, in addition to smaller amounts of MbA, the therapeutically inactive montbretin B (MbB) was the major product of this metabolic engineering effort. MbA and MbB differ in a single hydroxyl group of their acyl side chains, which are derived from caffeoyl-CoA and coumaroyl-CoA, respectively. Biosynthesis of both MbA and MbB also require coumaroyl-CoA for the formation of the myricetin core. Caffeoyl-CoA and coumaroyl-CoA are formed in the central phenylpropanoid pathway by acyl activating enzymes (AAEs) known as 4-coumaroyl-CoA ligases (4CLs). Here we investigated a small family of montbretia AAEs and 4CLs, and their possible contribution to montbretin biosynthesis. RESULTS: Transcriptome analysis for gene expression patterns related to montbretin biosynthesis identified eight different montbretia AAEs belonging to four different clades. Enzyme characterization identified 4CL activity for two clade IV members, Cc4CL1 and Cc4CL2, converting different hydroxycinnamic acids into the corresponding CoA thioesters. Both enzymes preferred coumaric acid over caffeic acid as a substrate in vitro. While expression of montbretia AAEs did not enhance MbA biosynthesis in N. benthamiana, we demonstrated that both Cc4CLs can be used to activate coumaric and caffeic acid towards flavanone biosynthesis in yeast (Saccharomyces cerevisiae). CONCLUSIONS: Montbretia expresses two functional 4CLs, but neither of them is specific for the formation of caffeoyl-CoA. Based on differential expression analysis and phylogeny Cc4CL1 is most likely involved in MbA biosynthesis, while Cc4CL2 may contribute to lignin biosynthesis. Both Cc4CLs can be used for flavanone production to support metabolic engineering of MbA in yeast.

Complete Biosynthesis of the Anti-Diabetic Plant Metabolite Montbretin A.

Diabetes and obesity are affecting human health worldwide. Their occurrence is increasing with lifestyle choices, globalization of food systems, and economic development. The specialized plant metabolite montbretin A (MbA) is being developed as an antidiabetes and antiobesity treatment due to its potent and specific inhibition of the human pancreatic α-amylase. MbA is a complex acylated flavonol glycoside formed in small amounts in montbretia (Crocosmia × crocosmiiflora) corms during the early summer. The spatial and temporal patterns of MbA accumulation limit its supply for drug development and application. We are exploring MbA biosynthesis to enable metabolic engineering of this rare and valuable compound. Genes and enzymes for the first four steps of MbA biosynthesis, starting from the flavonol precursor myricetin, have recently been identified. Here, we describe the gene discovery and functional characterization of the final two enzymes of MbA biosynthesis. The UDP-glycosyltransferases, CcUGT4 and CcUGT5, catalyze consecutive reactions in the formation of the disaccharide moiety at the 4'-hydroxy position of the MbA flavonol core. CcUGT4 is a flavonol glycoside 4'-O-xylosyltransferase that acts on the second to last intermediate (MbA-XR2) in the pathway. CcUGT5 is a flavonol glycoside 1,4-rhamnosyltransferase that converts the final intermediate (MbA-R2) to complete the MbA molecule. Both enzymes belong to the UGT family d-clade and are specific for flavonol glycosides and their respective sugar donors. This study concludes the discovery of the MbA biosynthetic pathway and provides the complete set of genes to engineer MbA biosynthesis. We demonstrate successful reconstruction of MbA biosynthesis in Nicotiana benthamiana.

Terpene Synthases and Terpene Variation in Cannabis sativa.

Cannabis (Cannabis sativa) resin is the foundation of a multibillion dollar medicinal and recreational plant bioproducts industry. Major components of the cannabis resin are the cannabinoids and terpenes. Variations of cannabis terpene profiles contribute much to the different flavor and fragrance phenotypes that affect consumer preferences. A major problem in the cannabis industry is the lack of proper metabolic characterization of many of the existing cultivars, combined with sometimes incorrect cultivar labeling. We characterized foliar terpene profiles of plants grown from 32 seed sources and found large variation both within and between sets of plants labeled as the same cultivar. We selected five plants representing different cultivars with contrasting terpene profiles for clonal propagation, floral metabolite profiling, and trichome-specific transcriptome sequencing. Sequence analysis of these five cultivars and the reference genome of cv Purple Kush revealed a total of 33 different cannabis terpene synthase (CsTPS) genes, as well as variations of the CsTPS gene family and differential expression of terpenoid and cannabinoid pathway genes between cultivars. Our annotation of the cv Purple Kush reference genome identified 19 complete CsTPS gene models, and tandem arrays of isoprenoid and cannabinoid biosynthetic genes. An updated phylogeny of the CsTPS gene family showed three cannabis-specific clades, including a clade of sesquiterpene synthases within the TPS-b subfamily that typically contains mostly monoterpene synthases. The CsTPSs described and functionally characterized here include 13 that had not been previously characterized and that collectively explain a diverse range of cannabis terpenes.

Biosynthesis of the anti-diabetic metabolite montbretin A: glucosylation of the central intermediate mini-MbA.

Type 2 diabetes (T2D) affects over 320 million people worldwide. Healthy lifestyles, improved drugs and effective nutraceuticals are different components of a response against the growing T2D epidemic. The specialized metabolite montbretin A (MbA) is being developed for treatment of T2D and obesity due to its unique pharmacological activity as a highly effective and selective inhibitor of the human pancreatic α-amylase. MbA is an acylated flavonol glycoside found in small amounts in montbretia (Crocosmia × crocosmiiflora) corms. MbA cannot be obtained in sufficient quantities for drug development from its natural source or by chemical synthesis. To overcome these limitations through metabolic engineering, we are investigating the genes and enzymes of MbA biosynthesis. We previously reported the first three steps of MbA biosynthesis from myricetin to myricetin 3-O-(6'-O-caffeoyl)-glucosyl rhamnoside (mini-MbA). Here, we describe the sequence of reactions from mini-MbA to MbA, and the discovery and characterization of the gene and enzyme responsible for the glucosylation of mini-MbA. The UDP-dependent glucosyltransferase CcUGT3 (UGT703E1) catalyzes the 1,2-glucosylation of mini-MbA to produce myricetin 3-O-(glucosyl-6'-O-caffeoyl)-glucosyl rhamnoside. Co-expression of CcUGT3 with genes for myricetin and mini-MbA biosynthesis in Nicotiana benthamiana validated its biological function and expanded the set of genes available for metabolic engineering of MbA.

Flavonol Biosynthesis Genes and Their Use in Engineering the Plant Antidiabetic Metabolite Montbretin A.

The plant metabolite montbretin A (MbA) and its precursor mini-MbA are potential new drugs for treating type 2 diabetes. These complex acylated flavonol glycosides only occur in small amounts in the corms of the ornamental plant montbretia (Crocosmia × crocosmiiflora). Our goal is to metabolically engineer Nicotiana benthamiana using montbretia genes to achieve increased production of mini-MbA and MbA. Two montbretia UDP-dependent glycosyltransferases (UGTs), CcUGT1 and CcUGT2, catalyze the formation of the first two pathway-specific intermediates in MbA biosynthesis, myricetin 3-O-rhamnoside and myricetin 3-O-glucosyl rhamnoside. In previous work, expression of these UGTs in N. benthamiana resulted in small amounts of kaempferol glycosides but not myricetin glycosides, suggesting that myricetin was limiting. Here, we investigated montbretia genes and enzymes of flavonol biosynthesis to enhance myricetin formation in N. benthamiana We characterized two flavanone hydroxylases, a flavonol synthase, a flavonoid 3'-hydroxylase (F3'H), and a flavonoid 3'5'-hydroxylase (F3'5'H). Montbretia flavonol synthase converted dihydromyricetin into myricetin. Unexpectedly, montbretia F3'5'H shared higher sequence relatedness with F3'Hs in the CYP75B subfamily of cytochromes P450 than with those with known F3'5'H activity. Transient expression of combinations of montbretia flavonol biosynthesis genes and a montbretia MYB transcription factor in N. benthamiana resulted in availability of myricetin for MbA biosynthesis. Transient coexpression of montbretia flavonol biosynthesis genes combined with CcUGT1 and CcUGT2 in N. benthamiana resulted in 2 mg g-1 fresh weight of the MbA pathway-specific compound myricetin 3-O-glucosyl rhamnoside. Additional expression of the montbretia acyltransferase CcAT1 led to detectable levels of mini-MbA in N. benthamiana.

Discovery of UDP-Glycosyltransferases and BAHD-Acyltransferases Involved in the Biosynthesis of the Antidiabetic Plant Metabolite Montbretin A.

Plant specialized metabolism serves as a rich resource of biologically active molecules for drug discovery. The acylated flavonol glycoside montbretin A (MbA) and its precursor myricetin 3-O-(6'-O-caffeoyl)-glucosyl rhamnoside (mini-MbA) are potent inhibitors of human pancreatic α-amylase and are being developed as drug candidates to treat type-2 diabetes. MbA occurs in corms of the ornamental plant montbretia (Crocosmia x crocosmiiflora), but a system for large-scale MbA production is currently unavailable. Biosynthesis of MbA from the flavonol myricetin and MbA accumulation occur during early stages of corm development. We established myricetin 3-O-rhamnoside (MR), myricetin 3-O-glucosyl rhamnoside (MRG), and mini-MbA as the first three intermediates of MbA biosynthesis. Contrasting the transcriptomes of young and old corms revealed differentially expressed UDP-sugar-dependent glycosyltransferases (UGTs) and BAHD-acyltransferases (BAHD-ATs). UGT77B2 and UGT709G2 catalyze the consecutive glycosylation of myricetin to produce MR and of MR to give MRG, respectively. In addition, two BAHD-ATs, CcAT1 and CcAT2, catalyze the acylation of MRG to complete the formation of mini-MbA. Transcript profiles of UGT77B2, UGT709G2, CcAT1, and CcAT2 during corm development matched the metabolite profile of MbA accumulation. Expression of these enzymes in wild tobacco (Nicotiana benthamiana) resulted in the formation of a surrogate mini-MbA, validating the potential for metabolic engineering of mini-MbA in a heterologous plant system.

A Conifer UDP-Sugar Dependent Glycosyltransferase Contributes to Acetophenone Metabolism and Defense against Insects.

Acetophenones are phenolic compounds involved in the resistance of white spruce (Picea glauca) against spruce budworm (Choristoneura fumiferiana), a major forest pest in North America. The acetophenones pungenol and piceol commonly accumulate in spruce foliage in the form of the corresponding glycosides, pungenin and picein. These glycosides appear to be inactive against the insect but can be cleaved by a spruce ß-glucosidase, PgßGLU-1, which releases the active aglycons. The reverse glycosylation reaction was hypothesized to involve a family 1 UDP-sugar dependent glycosyltransferase (UGT) to facilitate acetophenone accumulation in the plant. Metabolite and transcriptome profiling over a developmental time course of white spruce bud burst and shoot growth revealed two UGTs, PgUGT5 and PgUGT5b, that glycosylate pungenol. Recombinant PgUGT5b enzyme produced mostly pungenin, while PgUGT5 produced mostly isopungenin. Both UGTs also were active in vitro on select flavonoids. However, the context of transcript and metabolite accumulation did not support a biological role in flavonoid metabolism but correlated with the formation of pungenin in growing shoots. Transcript levels of PgUGT5b were higher than those of PgUGT5 in needles across different genotypes of white spruce. These results support a role of PgUGT5b in the biosynthesis of the glycosylated acetophenone pungenin in white spruce.

In vivo function of Pgßglu-1 in the release of acetophenones in white spruce.

Eastern spruce budworm (Choristoneura fumiferiana Clemens) (ESBW) is a major forest pest which feeds on young shoots of white spruce (Picea glauca) and can cause landscape level economic and ecological losses. Release of acetophenone metabolites, piceol and pungenol, from their corresponding glycosides, picein and pungenin, can confer natural resistance of spruce to ESBW. A beta-glucosidase gene, Pgßglu-1, was recently discovered and the encoded enzyme was characterized in vitro to function in the release of the defensive acetophenone aglycons. Here we describe overexpression of Pgßglu-1 in a white spruce genotype whose metabolome contains the glucosylated acetophenones, but no detectable amounts of the aglycons. Transgenic overexpression of Pgßglu-1 resulted in release of the acetophenone aglycons in planta. This work provides in vivo evidence for the function of Pgßglu-1.

Biosynthesis of Sandalwood Oil: Santalum album CYP76F cytochromes P450 produce santalols and bergamotol.

Sandalwood oil is one of the world's most highly prized essential oils, appearing in many high-end perfumes and fragrances. Extracted from the mature heartwood of several Santalum species, sandalwood oil is comprised mainly of sesquiterpene olefins and alcohols. Four sesquiterpenols, α-, ß-, and epi-ß-santalol and α-exo-bergamotol, make up approximately 90% of the oil of Santalum album. These compounds are the hydroxylated analogues of α-, ß-, and epi-ß-santalene and α-exo-bergamotene. By mining a transcriptome database of S. album for candidate cytochrome P450 genes, we cloned and characterized cDNAs encoding a small family of ten cytochrome P450-dependent monooxygenases annotated as SaCYP76F37v1, SaCYP76F37v2, SaCYP76F38v1, SaCYP76F38v2, SaCYP76F39v1, SaCYP76F39v2, SaCYP76F40, SaCYP76F41, SaCYP76F42, and SaCYP76F43. Nine of these genes were functionally characterized using in vitro assays and yeast in vivo assays to encode santalene/bergamotene oxidases and bergamotene oxidases. These results provide a foundation for production of sandalwood oil for the fragrance industry by means of metabolic engineering, as demonstrated with proof-of-concept formation of santalols and bergamotol in engineered yeast cells, simultaneously addressing conservation challenges by reducing pressure on supply of sandalwood from native forests.

Transcriptome resources and functional characterization of monoterpene synthases for two host species of the mountain pine beetle, lodgepole pine (Pinus contorta) and jack pine (Pinus banksiana).

BACKGROUND: The mountain pine beetle (MPB, Dendroctonus ponderosae) epidemic has affected lodgepole pine (Pinus contorta) across an area of more than 18 million hectares of pine forests in western Canada, and is a threat to the boreal jack pine (Pinus banksiana) forest. Defence of pines against MPB and associated fungal pathogens, as well as other pests, involves oleoresin monoterpenes, which are biosynthesized by families of terpene synthases (TPSs). Volatile monoterpenes also serve as host recognition cues for MPB and as precursors for MPB pheromones. The genes responsible for terpene biosynthesis in jack pine and lodgepole pine were previously unknown. RESULTS: We report the generation and quality assessment of assembled transcriptome resources for lodgepole pine and jack pine using Sanger, Roche 454, and Illumina sequencing technologies. Assemblies revealed transcripts for approximately 20,000 - 30,000 genes from each species and assembly analyses led to the identification of candidate full-length prenyl transferase, TPS, and P450 genes of oleoresin biosynthesis. We cloned and functionally characterized, via expression of recombinant proteins in E. coli, nine different jack pine and eight different lodgepole pine mono-TPSs. The newly identified lodgepole pine and jack pine mono-TPSs include (+)-α-pinene synthases, (-)-α-pinene synthases, (-)-ß-pinene synthases, (+)-3-carene synthases, and (-)-ß-phellandrene synthases from each of the two species. CONCLUSION: In the absence of genome sequences, transcriptome assemblies are important for defence gene discovery in lodgepole pine and jack pine, as demonstrated here for the terpenoid pathway genes. The product profiles of the functionally annotated mono-TPSs described here can account for the major monoterpene metabolites identified in lodgepole pine and jack pine.

Evolution of conifer diterpene synthases: diterpene resin acid biosynthesis in lodgepole pine and jack pine involves monofunctional and bifunctional diterpene synthases.

Diterpene resin acids (DRAs) are major components of pine (Pinus spp.) oleoresin. They play critical roles in conifer defense against insects and pathogens and as a renewable resource for industrial bioproducts. The core structures of DRAs are formed in secondary (i.e. specialized) metabolism via cycloisomerization of geranylgeranyl diphosphate (GGPP) by diterpene synthases (diTPSs). Previously described gymnosperm diTPSs of DRA biosynthesis are bifunctional enzymes that catalyze the initial bicyclization of GGPP followed by rearrangement of a (+)-copalyl diphosphate intermediate at two discrete class II and class I active sites. In contrast, similar diterpenes of gibberellin primary (i.e. general) metabolism are produced by the consecutive activity of two monofunctional class II and class I diTPSs. Using high-throughput transcriptome sequencing, we discovered 11 diTPS from jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta). Three of these were orthologous to known conifer bifunctional levopimaradiene/abietadiene synthases. Surprisingly, two sets of orthologous PbdiTPSs and PcdiTPSs were monofunctional class I enzymes that lacked functional class II active sites and converted (+)-copalyl diphosphate, but not GGPP, into isopimaradiene and pimaradiene as major products. Diterpene profiles and transcriptome sequences of lodgepole pine and jack pine are consistent with roles for these diTPSs in DRA biosynthesis. The monofunctional class I diTPSs of DRA biosynthesis form a new clade within the gymnosperm-specific TPS-d3 subfamily that evolved from bifunctional diTPS rather than monofunctional enzymes (TPS-c and TPS-e) of gibberellin metabolism. Homology modeling suggested alterations in the class I active site that may have contributed to their functional specialization relative to other conifer diTPSs.

Transcriptome mining, functional characterization, and phylogeny of a large terpene synthase gene family in spruce (Picea spp.).

BACKGROUND: In conifers, terpene synthases (TPSs) of the gymnosperm-specific TPS-d subfamily form a diverse array of mono-, sesqui-, and diterpenoid compounds, which are components of the oleoresin secretions and volatile emissions. These compounds contribute to defence against herbivores and pathogens and perhaps also protect against abiotic stress. RESULTS: The availability of extensive transcriptome resources in the form of expressed sequence tags (ESTs) and full-length cDNAs in several spruce (Picea) species allowed us to estimate that a conifer genome contains at least 69 unique and transcriptionally active TPS genes. This number is comparable to the number of TPSs found in any of the sequenced and well-annotated angiosperm genomes. We functionally characterized a total of 21 spruce TPSs: 12 from Sitka spruce (P. sitchensis), 5 from white spruce (P. glauca), and 4 from hybrid white spruce (P. glauca × P. engelmannii), which included 15 monoterpene synthases, 4 sesquiterpene synthases, and 2 diterpene synthases. CONCLUSIONS: The functional diversity of these characterized TPSs parallels the diversity of terpenoids found in the oleoresin and volatile emissions of Sitka spruce and provides a context for understanding this chemical diversity at the molecular and mechanistic levels. The comparative characterization of Sitka spruce and Norway spruce diterpene synthases revealed the natural occurrence of TPS sequence variants between closely related spruce species, confirming a previous prediction from site-directed mutagenesis and modelling.

An integrated genomic, proteomic and biochemical analysis of (+)-3-carene biosynthesis in Sitka spruce (Picea sitchensis) genotypes that are resistant or susceptible to white pine weevil.

Conifers are extremely long-lived plants that have evolved complex chemical defenses in the form of oleoresin terpenoids to resist attack from pathogens and herbivores. In these species, terpenoid diversity is determined by the size and composition of the terpene synthase (TPS) gene family and the single- and multi-product profiles of these enzymes. The monoterpene (+)-3-carene is associated with resistance of Sitka spruce (Picea sitchensis) to white pine weevil (Pissodes strobi). We used a combined genomic, proteomic and biochemical approach to analyze the (+)-3-carene phenotype in two contrasting Sitka spruce genotypes. Resistant trees produced significantly higher levels of (+)-3-carene than susceptible trees, in which only trace amounts were detected. Biosynthesis of (+)-3-carene is controlled, at the genome level, by a small family of closely related (+)-3-carene synthase (PsTPS-3car) genes (82-95% amino acid sequence identity). Transcript profiling identified one PsTPS-3car gene (PsTPS-3car1) that is expressed in both genotypes, one gene (PsTPS-3car2) that is expressed only in resistant trees, and one gene (PsTPS-3car3) that is expressed only in susceptible trees. The PsTPS-3car2 gene was not detected in genomic DNA of susceptible trees. Target-specific selected reaction monitoring confirmed this pattern of differential expression of members of the PsTPS-3car family at the proteome level. Kinetic characterization of the recombinant PsTPS-3car enzymes identified differences in the activities of PsTPS-3car2 and PsTPS-3car3 as a factor contributing to the different (+)-3-carene profiles of resistant and susceptible trees. In conclusion, variation of the (+)-3-carene phenotype is controlled by copy number variation of PsTPS-3car genes, variation of gene and protein expression, and variation in catalytic efficiencies.

Identification and functional characterization of monofunctional ent-copalyl diphosphate and ent-kaurene synthases in white spruce reveal different patterns for diterpene synthase evolution for primary and secondary metabolism in gymnosperms.

The biosynthesis of the tetracyclic diterpene ent-kaurene is a critical step in the general (primary) metabolism of gibberellin hormones. ent-Kaurene is formed by a two-step cyclization of geranylgeranyl diphosphate via the intermediate ent-copalyl diphosphate. In a lower land plant, the moss Physcomitrella patens, a single bifunctional diterpene synthase (diTPS) catalyzes both steps. In contrast, in angiosperms, the two consecutive cyclizations are catalyzed by two distinct monofunctional enzymes, ent-copalyl diphosphate synthase (CPS) and ent-kaurene synthase (KS). The enzyme, or enzymes, responsible for ent-kaurene biosynthesis in gymnosperms has been elusive. However, several bifunctional diTPS of specialized (secondary) metabolism have previously been characterized in gymnosperms, and all known diTPSs for resin acid biosynthesis in conifers are bifunctional. To further understand the evolution of ent-kaurene biosynthesis as well as the evolution of general and specialized diterpenoid metabolisms in gymnosperms, we set out to determine whether conifers use a single bifunctional diTPS or two monofunctional diTPSs in the ent-kaurene pathway. Using a combination of expressed sequence tag, full-length cDNA, genomic DNA, and targeted bacterial artificial chromosome sequencing, we identified two candidate CPS and KS genes from white spruce (Picea glauca) and their orthologs in Sitka spruce (Picea sitchensis). Functional characterization of the recombinant enzymes established that ent-kaurene biosynthesis in white spruce is catalyzed by two monofunctional diTPSs, PgCPS and PgKS. Comparative analysis of gene structures and enzyme functions highlights the molecular evolution of these diTPSs as conserved between gymnosperms and angiosperms. In contrast, diTPSs for specialized metabolism have evolved differently in angiosperms and gymnosperms.

Poplar defense against insects: genome analysis, full-length cDNA cloning, and transcriptome and protein analysis of the poplar Kunitz-type protease inhibitor family.

*Kunitz protease inhibitors (KPIs) feature prominently in poplar defense responses against insects. The increasing availability of genomics resources enabled a comprehensive analysis of the poplar (p)KPI family. *Using genome analysis, expressed sequence tag (EST) mining and full-length (FL)cDNA cloning we established an inventory and phylogeny of pKPIs. Microarray and real-time PCR analyses were used to profile pKPI gene expression following real or simulated insect attack. Proteomics of insect midgut content was used to monitor stability of pKPI protein. *We identified 31 pKPIs in the genome and validated gene models by EST mining and cloning of 41 unique FLcDNAs. Genome organization of the pKPI family, with six poplar-specific subfamilies, suggests that tandem duplications have played a major role in its expansion. pKPIs are expressed throughout the plant and many are strongly induced by insect attack, although insect-specific signals seem initially to suppress the tree pKPI response. We found substantial peptide coverage for a potentially intact pKPI protein in insect midgut after eating poplar leaves. *These results highlight the complexity of an important defense gene family in poplar with regard to gene family size, differential constitutive and insect-induced gene expression, and resilience of at least one pKPI protein to digestion by herbivores.

Dirigent proteins in conifer defense II: Extended gene discovery, phylogeny, and constitutive and stress-induced gene expression in spruce (Picea spp.).

Analysis of expressed sequence tags (ESTs) and full-length (FL)cDNAs from species of spruce (Picea spp.) revealed a family of 35 unique dirigent proteins (DIR) and DIR-like proteins. Phylogenetic analysis indicates the spruce DIR and DIR-like genes cluster into three distinct subfamilies, DIR-a, DIR-b/d, and DIR-f, of a larger plant DIR and DIR-like gene family. Gene-specific primers were designed for 31 unique spruce DIR family genes, and closely related isoforms, and used to evaluate patterns of constitutive expression, as well as responses to herbivory by stem-boring insects (i.e., white pine weevil, Pissodes strobi) in bark tissue and defoliating insects (i.e., western spruce budworm, Choristoneura occidentalis) in green apical shoots. Furthermore, meta-analysis of microarray gene expression data obtained from a series of independent experiments using the same 16.7K cDNA array platform identified several distinct expression clusters of the spruce DIR transcriptome closely matching phylogenetic clusters of sequence similarity. Members of the DIR-a family, which also contains functionally characterized DIR from other plant species, are most prominent for their induced response to feeding by weevils on Sitka spruce bark.

Aminocyclopropane carboxylic acid synthase is a regulated step in ethylene-dependent induced conifer defense. Full-length cDNA cloning of a multigene family, differential constitutive, and wound- and insect-induced expression, and cellular and subcellular localization in spruce and Douglas fir.

In conifer stems, formation of chemical defenses against insects or pathogens involves specialized anatomical structures of the phloem and xylem. Oleoresin terpenoids are formed in resin duct epithelial cells and phenolics accumulate in polyphenolic parenchyma cells. Ethylene signaling has been implicated in the induction of these chemical defenses. Recently, we reported the cloning of 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) from spruce (Picea spp.) and Douglas fir (Pseudotsuga menziesii). ACO protein was constitutively expressed in Douglas fir and only weakly induced upon wounding. We now cloned seven full-length and one near full-length cDNA representing four distinct 1-aminocyclopropane-1-carboxylic acid synthases (ACS; ACS1, ACS2, ACS3, and ACS4) from spruce and Douglas fir. Cloning of ACS has not previously been reported for any gymnosperm. Using gene-specific, quantitative real-time polymerase chain reaction, we measured constitutive expression for the four ACS genes and the single-copy ACO gene in various tissues of Sitka spruce (Picea sitchensis) and in white spruce (Picea glauca) somatic embryos. ACO and ACS4 were ubiquitously expressed at high levels; ACS1 was predominantly expressed in developing embryos and ACS2 and ACS3 were expressed only at very low levels. Insect attack or mechanical wounding caused strong induction of ACS2 and ACS3 in Sitka spruce bark, a moderate increase in ACO transcripts, but had no effect on ACS1 and ACS4. ACS protein was also strongly induced following mechanical wounding in Douglas fir and was highly abundant in resin duct epithelial cells and polyphenolic parenchyma cells. These results suggest that ACS, but not ACO, is a regulated step in ethylene-induced conifer defense.

Genomics of hybrid poplar (Populus trichocarpax deltoides) interacting with forest tent caterpillars (Malacosoma disstria): normalized and full-length cDNA libraries, expressed sequence tags, and a cDNA microarray for the study of insect-induced defences in poplar.

As part of a genomics strategy to characterize inducible defences against insect herbivory in poplar, we developed a comprehensive suite of functional genomics resources including cDNA libraries, expressed sequence tags (ESTs) and a cDNA microarray platform. These resources are designed to complement the existing poplar genome sequence and poplar (Populus spp.) ESTs by focusing on herbivore- and elicitor-treated tissues and incorporating normalization methods to capture rare transcripts. From a set of 15 standard, normalized or full-length cDNA libraries, we generated 139,007 3'- or 5'-end sequenced ESTs, representing more than one-third of the c. 385,000 publicly available Populus ESTs. Clustering and assembly of 107,519 3'-end ESTs resulted in 14,451 contigs and 20,560 singletons, altogether representing 35,011 putative unique transcripts, or potentially more than three-quarters of the predicted c. 45,000 genes in the poplar genome. Using this EST resource, we developed a cDNA microarray containing 15,496 unique genes, which was utilized to monitor gene expression in poplar leaves in response to herbivory by forest tent caterpillars (Malacosoma disstria). After 24 h of feeding, 1191 genes were classified as up-regulated, compared to only 537 down-regulated. Functional classification of this induced gene set revealed genes with roles in plant defence (e.g. endochitinases, Kunitz protease inhibitors), octadecanoid and ethylene signalling (e.g. lipoxygenase, allene oxide synthase, 1-aminocyclopropane-1-carboxylate oxidase), transport (e.g. ABC proteins, calreticulin), secondary metabolism [e.g. polyphenol oxidase, isoflavone reductase, (-)-germacrene D synthase] and transcriptional regulation [e.g. leucine-rich repeat transmembrane kinase, several transcription factor classes (zinc finger C3H type, AP2/EREBP, WRKY, bHLH)]. This study provides the first genome-scale approach to characterize insect-induced defences in a woody perennial providing a solid platform for functional investigation of plant-insect interactions in poplar.